ON THE MEASUREMENT OF KINETIC 
TEMPERATURE IN HIGH REDSHIFT GALACTIC 

HALOS 



j>^. Sergei A. Levshakov ^, Wilhelm H. Kegel ^ 

^ ' ^ Department of Theoretical Astrophysics, lojfe Institute, St. Petersburg, Russia. 

Institut fiir Theoretische Physik der Universitdt Frankfurt am Main, Frankfurt /Main, 

Germany. 



2 



' Abstract 

CN ■ We consider the deconvolution of the thermal and the turbulent profiles for a pair of 

ions with different atomic weights under the condition that the lines are formed in an 
■ absorbing medium with fluctuating density and velocity fields. Our method is based on 

On . the entropy-regularized x^-miniiiiization (ERM) procedure. From synthetic spectra we 

^ I found that the CII and Sill lines can be used to estimate the mean (density weighted) 

• kinetic temperature with a sufficiently high accuracy (~ 10%) in a rather wide range of 

the ionization parameter, —6.0 < logC/ < — 2.0. This is a generalization of a previous 
Q ■ result pi where we did not account for density fluctuations. 



1 Introduction 

One of the main source of information about the physical properties of the intergalactic gas is 
the study of metal-line absorption in quasar spectra. Metal lines in the Lyman Limit Systems 
(LLS), - i.e. in systems with Nuj < 2 x 10^'' cm~^ which are optically thin in the Lyman 
continuum, - are of particular interest in this regard because the Lyman continuum opacity 
does not affect the metal abundance measurements in this case [|l|]. Besides the LLSs often show 
carbon and silicon line absorption from different ionization states which allows estimating the 
ionization parameter U (the ratio of the number of photons with energies above one Rydberg 
to the number of atoms, U = n^/nn). 

The LLSs are usually assumed to arise in the outer regions of intervening galactic halos 
where the electron density is rather low, Ue ~ 10^^ — 10^^ cm^'^. For such tenuous gas the 
collisional ionization is not important in determining the ionization fractions of ions since typical 
kinetic temperatures for the gas showing absorption in CII - CIV lines and in Sill - SilV lines 
are of the order 10^ — 2 x 10^ K. It follows, that the outer parts of galactic halos are mainly 
photoionized and that the thermal and the ionization state of the gas may be specified by the 
ionization parameter U only, [|l|. Thus, for a given value of (or the specific radiation flux Jq 
at 1 Rydberg), the dispersion of U along the line of sight represents a varying gas density. 

In recent years, great efforts have been made towards the precise measurements of metal line 
profiles in QSO spectra obtained with high spectral resolution, FWHM ~ 5—7 km s~^. However, 



current theoretical models cannot give a complete prediction to match the observational data 

i- 

The observations often show complex structures of the line profiles. It is traditional to treat 
them using the standard Voigt fitting deconvolution procedure. This procedure is based on the 
assumption that the apparent fluctuations within the line profile are caused by density clumps 
('cloudlets') with different radial velocities. The non-thermal (turbulent) velocity field inside 
each cloudlet is accounted for in the so-called microturbulent approximation. It was shown in 
0, however, that the microturbulent analysis may produce unphysical kinetic temperatures. 

The more general mesoturbulent approach (see |Q and the references cited therein) is based 
on the assumption that the intensity fluctuations within the line profile arise mainly from the 
irregular Doppler shifts in the absorption coefficient caused by macroscopic large-scale, rather 
than thermal, motions. If the macroscopic velocity field has a correlation length not small as 
compared with the linear size of the absorbing region, then the radial velocity distribution may 
deviated significantly from the Gaussian model. 

In a previous paper we developed a method aimed at recovering the kinetic temperature 
from complex metal-line spectra assuming a homogeneous gas density and a random velocity 
field. Here we outline our first results for the case when both the density and the velocity fields 
are of random nature. 



2 The ERM procedure and results 

The entropy-regularized x^-minimization (ERM) procedure utilizes complex but similar ab- 
sorption line profiles of ions with different atomic weights to estimate the mean value of Tkm 
for the whole absorbing region. The similarity of the complex profiles of ions with different 
masses and ionization potentials stems from the equal ionization fraction for both of them. To 
illustrate this statement we used a model thoroughly described in |^ : an optically thin gas 
ionized by a typical QSO photoionizing spectrum p with the metallicity Z/Z0 = 0.1. 

The consequent steps of our computational experiment are shown in Fig. 1. Panel (a) 
presents the random velocity field u{x) = v{x)/aturb, where x is the space coordinate s along 
the line of sight in units of the linear size L of the absorbing region. The fluctuating density field 
n{x)/no is depicted in panel (b). Both the velocity and the gas density fields were calculated 
using the moving average method described in . To obtain the real gas density fluctuations, 
we used the log-normal distribution for the density contrast 6 = {n — no)/nQ with the rms 
value of cr^ = 1 and no = 2 x 10~^ cm~^. The rms value for the velocity field was assumed to 
be 20 km s^^, and the linear size of the region L = 5 kpc. The chosen specific radiation flux 
Jo = 5 X 10~^^ ergs cm~^ s~^ Hz~^ corresponds in our case to Uq = 1.3 x 10~^. In panel (c) we 
plot the equilibrium kinetic temperature T^ix) in units 10^ K as a function of the ionization 
parameter U{x) in accord with the numerical results |[l[]. The fluctuations of U{x) are shown 
in panel (d). For each species we calculated the density weighted mean temperature 

^ Jnion{s)T{s)ds 
/nion(s) ds 

The corresponding values of (T) we find for CII, Sill, CIV, and SilV are 15120 K, 15140 K, 
17120 K, and 16390 K, respectively. 

These random fields lead in turn to complex profiles of CIIA1334, SiIIA1260, CIVA1548, and 
SiIVA1394. All profiles have been convolved with a Gaussian spectrograph function having the 
width of 7 km s~^ (FWHM), and then a Gaussian noise of S/N = 75 has been added. The final 
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patterns are shown in panels (e), (f), (g), and (h), respectively, by dots with corresponding 
error bars. 

It is seen that the CII and Sill profiles are similar whereas the intensity fluctuations within 
the CIV profile differ from those for the SilV line. This different behavior of low and high 
ionized species formed in the same absorbing region is clearly illustrated in panel (i). The 
ratio of fractional ionization for the CII and Sill lines is almost constant over the whole range 
of ?7-values, but the corresponding ratio for the CIV and SilV lines is highly sensitive to the 
variation of U. 

Next we considered these spectra as 'observed' and analyzed them in two ways. At first we 
applied the standard Voigt profile fitting procedure. The results are shown in Figures le-lh, 
in where the individual components are indicated by tick marks and the numbers give the 
corresponding parameters (the first line gives the column density, the second the b parameter. 
In Fig. Ih the third line gives the derived temperature). The smooth curves show the resulting 
'theoretical' profiles. We see that the derived temperatures vary between 2800 K and 45900 K 
whereas the input temperature (Fig. Ic) varies only between 14360 K and 22090 K. Contrary, 
when we apply our new ERM procedure to the CII and Sill lines we find from Fig. Ij (shaded 
area) a temperature of 15200 ±200 K. This value is in good agreement with the density averaged 
temperatures derived before. 

This result reinforces our previous finding that the Voigt profile fitting procedure may lead 
to unphysical temperatures. On the other side it indicates that - at least for the cases studied 
so far - our new procedure gives a physically reasonable average value of the kinetic temperature. 

This is a report on work in progress. 
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Figure 1: An example of Monte-Carlo simulations, (a) One realization of the stochastic velocity 
field (cTturb = 20 km s~^), (6) the density field (no = 0.002 cm~^), (c) the kinetic temperature 
field (in units 10'^ K), and (d) the ionization parameter field {Uq = 1.3 x 10~^) versus the 
normalized space coordinate, (e, /, g, h) Velocity plots of the simulated spectra (dots and la 
error bars) and the Voigt profile fits to the CII, Sill, CIV, and SilV lines (smooth curves). 
The column density and velocity dispersion are displayed on top of each line. The third line 
in panel {h) shows the derived temperature for each subcomponent, {i) Fractional ionization 
ratios for the pair Cll/Sill and CIV/SilV lines versus the normalized ionization parameter, (j) 
The normalized cross-entropy (filled circles connected by smooth line), its curvature (dotted 
line) , and the kinetic temperature (dashed line) versus the normalized regularization parameter 
a (for more details see [Q). The shaded area restricts the region of optimal values of a and 
corresponding solutions for Tkin. 



